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We performed a series of 1341 full numerical simulations of high energy collision of black holes 
to search for the maximum recoil velocity after their merger. We consider equal mass binaries with 
opposite spins pointing along their orbital plane and perform a search of spin orientations, impact 
parameters, and initial linear momenta to find the maximum recoil for a given spin magnitude s. 
This spin sequence for s = 0.4, 0.7, 0.8, 0.85, 0.9 is then extrapolated to the extreme case, s = 1, to 
obtain an estimated maximum recoil velocity of 26,677 + 470 km/s, thus nearly 9% the speed of 


light. 
INTRODUCTION 


Ever since the discovery through full numerical sim- 
ulations [I] that the merger of binary black holes may 
lead to large (astrophysically speaking) gravitational re- 
coil velocities, a fascinating search for such events in na- 
ture takes place [2] [B]. Since the first modeling of large 
recoils [4], it was clear that the spins of the black holes 
played a crucial role in their merger remnant reaching up 
to several thousand km/s speeds. It was then found [5] 
a configuration that maximized the recoil nearing 5,000 
km/s. This configuration combined the opposite spins of 
that maximized asymmetry with the hangup effect [6] 
that maximized radiation. 


Those configurations assumed negligible eccentricities 
at the time of merger, when most of the asymmetric 
radiation takes place. While this is the most plausible 
astrophysical scenario, new gravitational waves observa- 
tions show the potential for large residual eccentricity 
in some events [7]. Here we will explore the extreme 
case of high energy collision of black holes, more in the 
realm of high-energy colliders and as tests of the radia- 
tion bounds theorems in General Relativity [8HI0]. This 
scenario was studied in [IT] to compute the maximum 
energy radiated by equal mass, nonspinning black holes 
in an ultrarelativistic headon collision. This first study 
was then followed up by the claim that the spin effects 
did not matter for these collisions in [I2]. 


In we revisited the headon scenario using our new 
initial data with low spurious initial radiation con- 
tent that allows for more accurate estimates of the max- 
imum energy radiated placing it at about 13%. Non 
headon high energy collisions have also been studied [15] 
and in notable detail in [16]. Some of the early reviews 
on the subject are and a more up-to-date one in 


Here we will focus our study on the computation of 
the maximum achievable gravitational recoil from a high 
energy collision of binary black holes and as we will see, 
the spin of the holes play a crucial role. 


NUMERICAL TECHNIQUES 


The full numerical simulations were performed us- 
ing the LAZEv code [20] implementation of the moving 
puncture approach [21]. We use the general relativis- 
tic BSSNOK formalism of evolutions systems [22}24]. 
The LAZEv code uses the Cactus [25] / CARPET [26] / 
EINSTEINTOOLKIT [27] [28] infrastructure. The CARPET 
mesh refinement driver provides a “moving boxes” style 
of mesh refinement. To compute the numerical (Bowen- 
York) initial data, we use the TWOPUNCTURES [29] code. 
We use AHFINDERDIRECT [30] to locate apparent hori- 
zons and measure the magnitude of the horizon spin 
Sy, using the isolated horizon algorithm as implemented 
in Ref. [BI]. We measure radiated energy, linear mo- 
mentum, and angular momentum, in terms of the ra- 
diative Weyl scalar 74, using the formulas provided in 
Refs. [32] [33]. As described in Ref. [34], we use the 
Teukolsky equation to analytically extrapolate expres- 
sions for ry, from a finite observer location to .4%* 

One can argue on asymmetry requirements that the 
maximum of the recoil should be obtained by equal mass, 
opposite spins configurations (hangup effect [6] here is 
replaced by the critical impact parameter separating di- 
rect merger to scattering). This configuration has been 
studied earlier in [35] [36] leading to a wide range of max- 
imum value estimates. We will revisit this problem with 
our new simulations to explicitly model the problem in 
terms of the initial velocity of the holes, v, impact pa- 
rameter, b, and spin, s = Sy/m%, (where my = mı, is 
the horizon mass of each hole). 


SIMULATIONS’ RESULTS 


Our simulations families consist of a choice of a spin 
magnitude, here s = 0.4,0.7,0.8,0.85,0.9, and for each 
of them an initial velocity (or momentum per mass), yv, 
and the (normalized by M) impact parameter, b, as mea- 
sured at the initial separation of the holes R = 50M 
(with M = mı + m2 the addition of the horizon masses 


V = (21799.58+/-331.54) sin( ọ + “5.654/-0.01 ) 
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FIG. 1. A run series versus ọ orientation of the spin for 


b = 2.38, s = 0.85, yv = 0.874. 
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FIG. 2. Maximum recoil velocity for different impact param- 
eters b. For all the high spins studied this peaks at b ~ 2.38. 


of the system). We then vary the orientation of the spins 
pointing on the orbital plane by an angle y with respect 
to the line initially joining the black holes. This allows 
us to model the leading y-dependence of the holes as a 
cos y [37]. In practice one needs about 5 to 7 simulations 
to fit this dependence and determine the amplitude of 
the curve leading to the maximum recoil configuration, 
as for instance displayed in Fig. 

This process is repeated for each impact parameter b 
to find the value bmaz that leads to the largest recoil ve- 
locity. Fig. |2| displays this search for each of the spin 
magnitudes we considered here. In practice, the bmaz 
corresponds closely to the critical value of the impact 
parameter b, separating the direct merger from the scat- 
tering of the holes. 

A similar analysis can be done by varying the initial 
velocity (or linear momentum) of the holes, yv, with y = 
(1 — v?)-'/?, the Lorentz factor. Those curves display 
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FIG. 3. Display of peak velocity vs. yv and spin for merging 
holes. 
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FIG. 4. Maximum recoil velocity for different initial momenta 
parameters yu. 


the same feature of maximizing the recoil velocity for 
values about the critical momentum separating the direct 
merger from scattering of the holes as displayed in Fig. 
The return of the curves represent an overshot of the 
impact parameter compensated by the lowering of the 
initial velocity to warrant merger. 

The explicit dependence on both parameters, b and yu 
is displayed in Fig. [4] as a heat map for each of the spin 
magnitudes s = 0.7, 0.8, 0.85, 0.9. 

The final results of the maximum recoil velocities for 
each s is summarized in Table [I] and in Fig. |5| where we 
display the error bars of each point and fit a quadratic 
dependence with s to extrapolate to the ultimate recoil 
velocity finding 26,677+470 km/s for the extremely spin- 
ning, s = 1, binary black holes case. 

Note also that in Table |I| we provide the number of 
runs (simulations) we performed in the three-dimensional 
search (y, b, yv) of the maximum recoil for a given spin s. 
In the case of s = 0.4, to simplify the search, and on the 
light of the previous results with the higher spin cases 
(as seen in Figs. 2]and|4) we assumed b = 2.38, hence the 


TABLE I. All simulations have equal mass q = mi/mz2 = 1, 
and placed at 21,2/M = +25 


+s aad (VW) maz Vmax [km/s] # runs 
0.40 2.38 1.20 11,590 + 50 64 
0.70 2.38 1.10 19,800 + 260 225 
0.80 2.38 1.10 21,900 + 210 464 
0.85 2.38 1.10 23,000 + 500 297 
0.90 2.38 1.09 24,100 + 590 291 


peak V [km/s] 


— V(s=1) = 26677.06 +/- 469.54 km/s 
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FIG. 5. Maximum recoil velocity versus spin of the holes in 
the q = 1 opposed in orbital plane configuration 


lower number of simulations needed. 

As a further check of our numerical accuracy, we have 
recalculated the cases of spin 0.8 with the extra refine- 
ment level and increased grid sizes we are using for the 
spin 0.9 runs. We then recalculated the series that gives a 
maximum value for spin 0.8 of 21802 + 191 km/s. Com- 
pared to the original grid computation of 21903 + 213 
km/s, leads to a difference of 101 km/s or 0.46%. 


CONCLUSIONS 


We have been able to provide an accurate estimate of 
the ultimate recoil, product of the high energy collision 
of black holes. In order to perform the four dimensional 
search (momentum yv, impact parameter b, spin orienta- 
tion y and magnitude s) we performed 1341 simulations 
in search for the critical be marginally leading to merger 
and the corresponding value of Ymax that maximized the 
recoil all as a function of y. Extrapolation to maximum 
spins have led us to estimate the value of 26,677 + 470 
km/s for the ultimate recoil, placing a bound below 10% 
of the speed of light. 

In Fig. [6] we display the spectrum of radiated energy 
adding the leading (¢ = 2,m = 0,+2)-modes for one 
of the peak recoil cases (b = 2.38, s = 0.85, yv = 1.1) 
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FIG. 6. The spectrum of the (£ = 2-modes) energy radiated 
dEv=2/dw by a representative set of simulations (with b = 
2.38, s = 0.85, yv = 1.1) for different orientation angles, p of 
the spin. 


for different orientation angles, y of the spin. We observe 
that at low frequencies, that correspond to the initial and 
“bremsstrahlung-like” radiation of the holes approaching 
each other, the spin orientation does not display notable 
differences, while at higher frequencies, corresponding to 
the holes merger, we observe a strong dependence on the 
spin orientations. 

We thus note the crucial relevance of the spins of the 
holes in the determination of the high energy collision 
kicks. In a follow up study we will determine in detail 
the role of the spins in the maximum radiated energy and 
angular momentum [38]. 
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